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Abstract 

In this paper we investigate spherically symmetric vacuum solutions of f(R) gravity in 
a higher dimensional spacetime. With this objective we construct a system of non-linear 
differential equations, whose solutions depend on the explicit form assumed for the function 
F(R) = j ■ We explicit show that for specific classes of this function exact solutions from 
the field equations are obtained; also we find approximated results for the metric tensor for 
more general cases admitting F(R) close to the unity. 
PACS numbers: 04.50. + h, 04.20. - q 
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oo ■ 1 Introduction 

O . 

One of most fascinating problem in modern Cosmology is to explain satisfactorily the accelerated 
expansion of the Universe pp. In this way three main candidates have been considered in 
literature: The presence of a cosmological constant A, the existence of a dark energy and the 
modified theories of gravity. The cosmological constant variant brings naturally the cosmological 
problem and the coincidence problem as well. Requiring that A presents a tiny value to cause the 
present acceleration demands to an extreme fine-tuning. In the second scenario, it is postulated 
the existence of a fluid with equation of state P ks —p, being p and P the energy density and 
pressure, respectively, which comes to dominate late in the matter era. The third alternative 
opens the possibility to provide new contributions to the equation of state compatible with an 
accelerated expansion of the Universe. The main idea of modified gravity theories resides in 
the generalization of the Einstein-Hilbert action in large scale, admitting a scalar non-linear 
function of the Ricci scalar curvature, /(J2)0 The modified gravity also allows a unification 
of the early-time inflation [3] and late-time cosmic speed-up 0J [5] . □ These models seem also 
relevant to explain the hierarchy problem and unification of Grand Unified Theories, GUTs, 
with gravity [8j. Such theories avoid the Ostrogradski's instability [9] that can otherwise prove 
to be problematic for general higher derivatives theories |10j . 

In the context of modified theories of gravity, static spherically symmetric solutions of the 
field equations have been analysed in vacuum sector [11] as well as in the presence of perfect fluid 
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The idea of modifying Einstein-Hilbert action has already been proposed many years ago [2], as a possible 
approach to construct a renormalized theory of gravity. 

Specifically in [6] [7], the authors proposed realistic models which support the inflationary epoch with the 
present cosmic acceleration. 
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source [12) . respectively; in the later, the authors show that a given distribution of matter is not 
capable to uniquely determine f(R)- In |13j , spherically symmetric solutions of f(R) theories 
of gravity have been analysed by using the Noether symmetry. Moreover, static cylindrically 
symmetric vacuum solution has been analysed in [14] . 

Another attractive model to describe our Universe is the braneworld scenario. By this 
model our world is represented by a four-dimensional sub-manifold, a three-brane, embedded in 
a higher dimensional spacetime [151 116] . In fact the conjecture that our Universe may have more 
than four dimensions was proposed by Kaluza [T7] with the objective to unify gauge theories 
with gravitation in a geometric formalism. In this context extra coordinates are compactified 
to circles of small radius that would be observable at high-energy scale comparable to the 
Planck one. Many of high energy theories of fundamental physics are formulated in higher- 
dimensional spacetime. In supergravity and superstring the idea of extra dimensions has been 
extensively used. Also in higher-dimensional spacetime context, of particular interest are the 
models introduced by Randall and Sundrum [181 EH] • I n these models it is assumed that all 
matter fields are confined on the brane and gravity only propagates in the five dimensional bulk. 
In the RSI model, the hierarchy problem between the Planck scale and the electroweak one is 
solved if the distance between the two branes is about 37 times the AdS radius. 

In this paper we join both promising ideas about the Universe. In this sense we shall analyse 
the modified theories of gravity in a higher dimensional spacetime under classical point of view. 
Specifically, as the first starting point, we shall analyse vacuum solutions admitting spherically 
symmetric Ansatz for the metric tensor. In this sense our approach is a higher dimensional 
extension of that one presented in [11]. We shall adopt radial functions for the function F(R) 
and calculate the corresponding results for the metric tensor. 

This paper is organized as follows. In Section [2] we present for a (1 + d)— dimensional 
spacetime, spherically symmetric fields equations for general modified theories of gravity. We 
also present some relevant properties related with the geometry of the corresponding spacetime. 
In Section [3] we provide explicitly vacuum solutions of the field equations for a large class of radial 
function F(R). We also analyse the specific case with d = 2. We show that differently from 
Einstein equation, non-trivial vacuum solutions are obtained for the latter. Moreover; comparing 
this case with large value of d ones, we show that there exist a substantial modification on the 
field equations, which justify a separated analysis. Finally we summarize our most important 
results in Section 0[ The Appendix [A] contains some technical details useful in Section [3l In this 
paper we use signature +2, and the definitions: Rp^s = ^7^/35 ~~ •••> R&P = Ra-y/3- We also use 
units with G = c = 1. 

2 Field equations in higher-dimensional spherically symmetric 
spacetime 

In this section we shall analyse the modified gravity action in a higher-dimensional spacetime. 
Specifically static spherically symmetric metric tensor will be considered. We shall present some 
geometric properties associated with this corresponding spacetime, and provide the general fields 
equations obeyed by the components of the metric tensor. 

Let us consider first the action below for modified theories of gravity: 



where R is the Ricci scalar, k = 8ttG, and S m represents the action associated with matter 
fields. By using the metric formalism, the field equations becomes: 
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in which T£ u is the geometric energy-momentum tensor, namely 



= j^u, (f(R) ~ F(R)R) + V a V?F(R) (g a , gf3u - g^g a A (3) 

with F(R) = 

The standard minimally coupled energy- momentum tensor T™, derived from the matter 
action, is related to T™ by 

T™ = f™/F(R) . (4) 
In absence of matter fields the field equations reads 

F(R)R^ - \f(R)g^u - V„V V F{R) + g^UF(R) = . (5) 
Taking the trace of the above equation we get 

F(R)R - ^^-f(R) + dOF(R) = , (6) 

which express a further degree of freedom that arises in the modified theory, namely the scalar 
curvature one. Through this equation it is possible to express f(R) in term of its derivatives, as 
follows 

f(R) = ^-(F(R)R + dOF(R)) . (7) 
Substituting the above expression into ([5]) we obtain 

F(R)R^ - V^V„F(R) = [F(R)R - DF(R)] . (8) 

From this expression we can see that the combination below 

A = F(R)R^ - Vj^HR) ^ 

with fixed indices, is independent of the corresponding index. 



2.1 Spherically symmetric metric tensor 

Adopting the coordinate system = (t, r, 0%, O2, #<2-2> </>)} with d > 3 and the coordinates 
defined in the intervals t S (— 00,00), B% £ [0,tt] for i = 1,2. ..a! — 2, <j) € [0, 2tt] and r > 0, the 
static metric tensor associated with the most general spherically symmetric (1 + d)— dimensional 
spacetime, can be represented by 

5oo = -u(r) 

gu = V (r) 
922 = r 2 

gjj = r 2 sin 2 6\ sin 2 62... sin 2 Oj-2 , (10) 

for 3 < j < d, and g^ v = for /j, 7^ v. For this metric tensor the non-vanishing components of 
the Ricci tensor read: 

{u') 2 u'v' u" (d-lVl 

- — - — I — , (11) 

2u 2 v 2uv 2 uv uvr 



R 
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u'v' {d-l)v' u" {u'f 



+ - — ^ + 



(12) 



2uv 2 rv 2 uv 2u 2 v 
J u' (d-2) (d-2) 

^2 - ^—o ~ 77— IT— H "5 ■ l 1 ^) 



The spherical symmetry requires that 

-^2 = ^3 = ^4 = ••• = R d i 

which implies that the scalar curvature will be given by 

_ 2 f d(d-3) + 2 a[d(d-3) + 2] (d - 1) . . 

r 2 u \ 2 2 + 2 V u u 

In' 2 /V iA 1 2 u" 



V - + - +-r^— ^ . (15) 

In all the above equations the primes corresponds derivative of the function with respect to the 
radial coordinate. 

Defining Y(r) = u(r)v(r) and by using the identity = A v for all /x and is, we can construct 
the two linearly independent homogeneous differential equations below: 

Y' Y' 
2rF"-r —F'-(d-l) — F = Q, (16) 



and 



F' F 1 Y 1 

- 4(d - 2)u + AY(d - 2) - 4ru— + 2rV— - rV— 

F FY 

+ 2r 2 u" + 2ru'(d - 3) + 2ur— = . (17) 
These equation will be the basis for out future analysis. 

3 Vacuum solutions 

In this section we shall analyse possible vacuum solutions of the field equations with radial 
symmetry. As we shall see, these solutions strongly dependent on the explicit form adopted for 
F(R). 

3.1 Spacetimes with constant scalar curvature 

As first part of our investigation, let us consider spacetimes with constant scalar curvature. For 
this case, the field equations (fl~6j) and (JTTj) can be written as 

uv' + vu = (18) 

and 

1 A ' 2\u v J \2 u J 2 u 2u y ' 

whose solutions are shown to be given by 



v ( r ) = ~~~7~~\ wrtn u(r) = co + + C2r 2 , (20) 
u(r) r a ~ z 
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where cq, C\ and c 2 are integration constants. In order to maintain unchanged the signature 
we impose that cq > 0. By using (|15p . we can find that the corresponding scalar curvature is 
R = —d(d + 1)02/ cq. It is possible to rescale the time coordinate so that we can always choose 
Co = 1. A Lagrangian in which f(R) = R+A gives rise to Schwarzschild-de Sitter (SdS) solution, 
that means a Schwarzschild spacetime in the presence of a cosmological constant. 
For a d— dimensional space, the Newtonian potential obeys the equation 

V 2 $ = n d M5 (d) (7) , (21) 

being M the mass of the particle placed at origin. However, in the case where d > 3 it is known 



tha" 
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V '^2 =-(d-2)« w «W(7) , (22) 



where = jfrar • So, by using the above relations this potential is written as 

Ml 

* - ■ (23) 

Now coming back to the SdS-type solution, we have the forms: 

1 2M 1 Ar 2 

v(r) = , -. with uir) = 1 — — , „ — , (24) 

v ' u(r) y ' (d-2)r d ~ 2 d(d - 2) ' y ' 

whose corresponding scalar curvature is 

(25) 

This solution is achieved in (f20|) by choosing c\ = —2M/(d — 2) and ci = —2A/d(d — 1). 
Because we are deal with spherically symmetric solutions, the idea of horizon associated with 
the corresponding black hole, takes place. These regions are found by imposing the condition 
u(r) = 0; however for a higher dimensional spacetime this condition leads us to high order 
algebraic equations whose solutions, in general, are not provided by any analytical procedure. 



3.2 Solutions with w(r)v(r)=constant 

Another class of solutions can be obtained by imposing Y(r) = u{r)v{r) = Yq =const. In this 
case a solution of (116j) is F(r) = Ar Substituting this expression into (I17p . we obtain an 

exact solution for u(r) given bjj^l 

v U A /2B 2 Y A \ {-l) d /A\ 3 - d 2 

a (~l) n 
B ^ nr n ~ 2 

n=3 

where Cj are integrations constants. In order to have a SdS-type solutions for the non-trivial 
gravity theory, the Newtonian potential- type term requires that c\ = — ^_^^. 2 ) (5) ' ^ e 
can define the independent term Yq + \ = 1, which leads Yq ^ 1. In this case, we shall have 
a correction in the Newtonian potential due the term Ac\/2B 2 7^ 0. It is natural to choose 

3 The d = 2 case, will be analysed separately in subsection 13.41 

4 In this paper we shall use F(r) to represent the radial function obtained from F(R(r)). 
5 The details of this calculation are in Appendix lAl 



(AX'-" ,,A' /2B'Y d , V / B\ 
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c\ = — (_tyj|^_2) i^) d 3 m or der the form of Newtonian potential be achieved. In this case we 
have: 



Yn 



and 



.4 



d- 1 
d-2 

(— l) d \ 
dM 



where we have set B = 1. This leads to the following form for jP(r) 



F(r) = 1 + 



dM 



r . 



(27) 
(28) 

(29) 



Finally, (|26|) will be given by 

2M 



u(r) = 1 



(d - 2)r d ~ 2 d(d - 1) 



2A „ 2dM t-k (-l) n + 1 - d 
r 2 + ^J2 



d — 2 < — ' nr 

n=3 



n-2 



(-1) 

dM 



d — n 

d\ d-2 



(30) 



Since Y(r) = u(r)v(r), through (j27[) we can straightforwardly determine v(r). The corresponding 
scalar curvature is given by: 



D d-2 2A(d-2)(d+l) ^M«^ 

R=^^ + — - — — — + 2^ — ^ — 



(d-lf 



(31) 



n=3 



which depends on the cosmological constant and the mass parameter. Moreover, one can obtain 
f(R) by using (fT6|) . So, the exact form, in terms of the radial coordinate, r, for f(r) will be 
given by 



f(r) 



(d+1) 



d-2 2A(d-2)(d+l) ^(-l) d V d_2) (d-2)(d + l) 



dM 



(d-l)r d ~ lV 7 ^ 



Md-2 Cn(J Md-ip n4 M*-*s ntd 



+ 



r n— 1 



+ 



„n— 1 



n=3 



(32) 



Where we have defined 



2d^2 
n(d- 1 



.(_l)3d-2(n+i) |[ 2 (n - 1) - d] (d - 1) - (n - 2)(n - 1)} 



71 — 3 

2d"rf 3 2 3d 2 -d(2n+7)+4(n+l) 

Pn,d = — tt (-1) 3=5 {[2( re " 1) " d] (d - 1) - (n - 2)(n - 1)} 



n(d- 1) 



and 



d-|-n — 5 

2d d " 2 , 2(2d-n-l) 

sn,d=— — tt(-1) d " 2 (d-n + 1). 



(33) 
(34) 

n(<t-iy - - (35) 

It is important to observe that the coefficients p n <2, s n,d an d c„ <j will be present only when d > 4, 
so p n ,3 = s nt 3 = c n ^ = 0, Vn. The expressions ([3T]) and ([32]) reproduce previous results when 
we substitute d = 3. However, in principle, for d > 3, it is not a trivial matter to provide an 
explicit expression for f(R) by using the above results. So, as an application of this formalism, 
let us consider d = 4. For this case, we have the following exact results as shown below: 



(36) 



R _ 2 20A 8VM 
r 2 9 9r 3 
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and 

4 8A 2 8\/M 

f{r) = ^- + T + 7m-^- (37) 

It is possible, by using (|36|) to express the radial coordinate, r, in terms of the scalar curvature, 
R; this provides a third order algebraic equation, which by its turn possesses one real root. 
Because the result found for this root is a very long one, the corresponding expression for f(R) 
does not clarify the physical content of the system under investigation, in this way we shall only 
analyse the above expressions in two asymptotic regimes; they are: (i) for r 3> yM and (ii) 
r <C Let us start with the first approximation. For this case the leading term in (136p 

reads: 

2 20A 

R ~ -o + -s- ' ( 38 ) 
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consequently 



2 /9i?-20A 8A , , 



On the other hand, the second approximation leads: 



8VM 20A 

«---^ + — m 

and 

(41) 

Finally we would like to say that, in principle, similar procedures as we have done before can 
also be applied for higher-dimensional spacetimes. 

3.3 Different Ansdtze for F(r) 

Although the most interesting models for modified gravity theories are those that reproduce 
Einstein theory of gravity for small disance, in this subsection we shall consider specific Ansdtze 
for F(r) that allow us to obtain exact solutions for the metric tensor. 

3.3.1 Solutions with F(r) = F Q r 

Considering F(r) = Fqt, Eq. (|17|) provides an exact solution for u(r). It is: 

u(r) = Y ^ + -^ + c^. (42) 

It would be interesting if we included in this analysis only purely de Sitter-type solutions. 
Pursuing this objective we may set c\ = 0, Yb = (rf_2) anc ^ ° 2 = ~ d(d-i) • This provides 

9A 

»M = 1 - w^rf ■ w 

The scalar curvature associated with this solution is 

(d+1) 2 + d(d-3) 

R -jd^T) K + ^2 . ( 44 ) 



3 For d = 4 the geometric mass has unity of length square. 
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which leads, by means of the same procedure adopted in the previous section, to the following 
exact form for f(R) 



f(R) = ± 2i W-3) / (2 + d(d-3))(d-l) 



(d-1) 2 V 2R(d - 1) - 2A(d + 1) 



[d(d - 3) + 2d(d - 2) + 2] F / 2(d-l)fi-2(d+l)A 

d+i y (2+d(d-3))(d-i) ■ ( j 

The corresponding result matches with that one obtained by T. Multamaki and I. Vilja in pT] 
for d = 3. 

3.3.2 Solutions with F(r) = F r n 

Another Ansatz that we can investigate is that one in which F(r) = For n , for h / 1, whose 
corresponding solution for u(r) is given by 

u(r) = u r m , (46) 

where m = 2 "f",~V ; furthermore we have 

Y (r) - (d 2 + d(2n-3) + 2(n-lf)(d-(n-ir) 
{) ~ (d-2)(n + d-l) 2 ' 1 j 

With the help of the equations ([4*6]) and (|i7j) we can compute the scalar curvature. It is: 

^ 2)(„ 2) 1 

n(n-2) - (d- 1) r 2 ' v ; 

from which we can obtain the corresponding f(R) given by 

M ; (2-n) \n(n-2) - (d-l) J V ; 

Where the constant h n ^ is defined by 

_ 2d 2 (n - 1) + 6n(n - 1) + d(d 2 - 1) - 2(n 3 - 1) 
n ' d = (d + l)(d(d - 3)) + 2n(d + n) - 2(2n - 1) ' (50) 

However, in an empty space this constant has no physical relevance, because it can be dropped 
out when the action is extremized. 

3.3.3 Solutions with F(r) = 1 + Ar n with Ar n < 1 

As we have already mentioned, in general it is a difficult matter to obtain an explicit form 
for f(R); the reason is because some expressions found for the scalar curvature do not allow to 
express the radial coordinate, r, algebraically in terms of R. This is what happens in the present 
analysis when we use F(r) = 1 + Ar n . A possible way to circumvent this problem, is to adopt an 
approximated procedure. In this direction we shall consider that the system under consideration 
is a small correction on General Relativity Theory (GR); then we assume F(r) = 1 + Ar n with 
Ar n < 1. In this case, we were able to find an approximated solution for u(r): 

,\ i 2M 1 2A 2 2Ar n ( Md \ 
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The expression for Y(r) can be obtained from F(r) by using (|16p . In this case, taking into 
account the adopted approximation, we have: 

Y(r) = 1 + 2 ^2i Arn + °(( Arn ) 2 ) ■ ( 52 ) 
consequently, the scalar curvature can also be found approximately by 



R 



d 



In the above expression we have defined = 2dn ^^ ^ , B^ d = _ 4A ( d (^^-|+ n +1 ) anc j 

= 4A M(3rf-^-4>-rf+i] ^ ^hese resu lts are valid for d ^ 1, 2. 

3.4 Special case d = 2 

Although the main objective of this paper is to investigate solutions of the modified theories of 
gravity in a higher-dimensional spacetime, in this subsection we shall analyse the special case 
of planar f(R) gravity with azimuthal symmetric spaces. The reason for this analysis resides in 
the fact that, in three-dimensions, the Einstein tensor, in a source-free space, is identically zero 
and only trivial solutions are possible. This fact happens because the Riemann tensor and the 
Ricci one have the same number of components. In fact, these tensors are related by [20\ I21j : 

Kb = ( R l ~ k^) • (54) 



However, we would like to emphasize that in the modified theories of gravity, non-trivial vacuum 
solutions of the field equations arise. Considering d = 2, the equations (fT6j) and (fT7|) can be 
written as 

2rF" - r—F' - —F = (55) 

and 

// 1 ,Y' , F' F' Y' , 
ru — -ru — + u r— — 2u— + u— — u = . (56) 

By solving algebraically (|55|) for Y' jY and substituting into ([56]) . we obtain a differential equa- 
tion relating F and u, namely 

\rF' + F F t) \rF> + F rF J ' K ; 

whose general solution u(r), for an arbitrary F(r), is given by 

r r 2 F"(r)F(r)-r 2 F'(r) 2 +F(r) 2 , 

u(r) = C2r 2 + c\r 2 J ^ dr . (58) 

3.4.1 Solutions with constant curvature 

The field equations for constant curvature in the specific case of d = 2, can be directly obtained 
from the general equations (fT8|) and (fTUI) : 

m; + vu = (59) 
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which provides the following solutions 

v ( r ) = ~~~7~~~\ an d n ( r ) = c o + c\r 2 . (61) 
u(r) 

The scalar curvature associated with this metric, obtained from (|15p . is R = — Qc\/cq. The 
scalar curvature, in (2+1) dimensions, assuming the existence of a cosmological constant A, is 
given by R = 6A. In this case, we should choose c\ = —A. Furthermore, it is convenient, for 
reasons previously presented the choice of co = 1. Prom this solution, we can notice that there 
is no Schwarzschild-type term in this (2+1) dimensional framework, for constant curvature. It 
was already expected because the different form of the Newtonian potential. 



3.4.2 Solutions with u(r)v{r) =constant 

Also for this case it is easy to verify by using (|55p that Y{r) = Yq, being Yq a constant, 
consequently F(r) = Ar + B. This form provides the following solution: 

u (r) = c 2 r 2 + ci [B{B-2Ar) + 2A 2 r 2 In (A + B/r)) . (62) 

Apart from the factor r 2 , the logarithmic term can be related with the Newtonian potential in 
two spatial dimensions, since in this case it has a logarithmic form. However, the details about 
the Newtonian potential is not our main goal, so it is interesting for our present purpose to 
find a way to avoid this logarithmic term. A natural way to remove it would be to consider 
A = 0, but this is a trivial case which does not bring new feature to our analysis. An alternative 
procedure that we are going to adopt is to consider this analysis as an approximated one. This 
approach will be taken into account in the next section. 



3.4.3 An approximated solution 

Here we shall consider the form F(r) = 1 + Ar; however, implicitly assuming that Ar < 1. This 
corresponds to a small perturbation on GR. In this approximation, all the quadratic terms in 
Ar will be neglected. So, the solution given by (|62p will be written as 



u(r) C2r 2 + ci — 2Ac±r . (63) 

As we have already seen, one can always set Yq = 1, consequently v(r) = l/u(r). We also 
consider c\ = 1 and admitting the existence of a cosmological constant, it is natural to choose 
the integration constant c 2 = —A. So we obtain: 

u{r) rs 1 - Ar 2 - 2Ar . (64) 

The approximated result for the scalar curvature obtained by substituting (|62p into (|15p is: 

4A 

R « 6A H + 12A 2 ln(r) + 1(L4 2 - 2AA 3 r . (65) 

r 

The corresponding form of f(R), in terms of the radial coordinate, r, is: 

, 4A A n, A . , 2(L4 2 28A 3 r 8AAr , . 

f(r) « 4A + — + 8^ 2 (1 + Ar) ln(r) + — — + — — . (66) 

T O O O 
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3.5 Solutions with F(r) = F r r 



Another Ansatz for F(r) that can also be examined is by considering F(r) = For n , with 

2ra(ra-l) 

For this case, equation (|55p provides the solution Y(r) = Y$r n+1 , being Yq a constant. 
Substituting this result into (f56|) . we obtain for u(r) the result below: 



u{r) = c\r 2 + c<ir n+1 . 



(67) 



Considering, as in the previous analysis, the coefficient of the term squared in the radial distance 
associated with a " cosmological constant"', we may take c\ = —A. Finally for this case the scalar 
curvature and the corresponding function f(R) read: 



R 



f(R) 



2{2n 2 - 5n - 3) A 



2n(n-l) 
f n + 1 



Y (n + 1) 



2F (n-l) 



n 



2(-2n 2 + 5n + 3)A 
Y (n + 1) 



n + 1 

2(n-l) 



^2(n-l) , 



(68) 



3.6 The Chern-Simons action 



In three dimensional spacetime, it is possible to consider the Chern-Simons topological mass 
term to analyse the dynamics associated with gravitational field [20|, 121] . In the context of f(R) 
theories of gravity, the most general action in (2 + 1)— dimensions has the following form: 



s= l.J d3xV =- sm+ ± Sc 



where 



(69) 



(70) 



being fi an arbitrary constant with dimension of mass, and lu^ the spin conection. The variation 
of S with respect to the metric tensor provides the field equation, 



Gov C uv 



(71) 



being T% v the geometric energy-momentum tensor, and C av the Cotton- York tensor defined by 

^ = V=SMV Q (<-W • (72) 

For the azimuthal symmetric space under consideration, the only non-vanishing component of 
the Cotton- York tensor is: 



i W 



02 



-U 



[6u"'Y 2 r 2 + 8u"Y 2 r - 8u'Y 2 + 4urY' 2 + 2uYY' + 6u'Y' 2 r 2 



8Y 4 r 

2urYY" - 9u"YY'r 2 - 8u'YY'r - 3u'YY"r 2 ] 



(73) 



which must vanish in the case of conformally flat space. We can verify that de Sitter type 
solution (|6T1) satisfies the condition C02 = 0, consequently is a solution of the field equation ([TTj) 
with constant curvature. 
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4 Conclusions 



In this paper we have analysed spherically symmetric vacuum solutions of f(R) theories of 
gravity in (1 + d)— dimensional spacetime. In order to achieve this objective, we have derived, 
by using the metric formalism, the field equations, which are expressed in terms of a general 
radial function F(r). Considering different models, i.e., different Ansatze for this function, 
exact and approximated solutions have been found. All the exact cases considered generalize 
previous results obtained by Multamaki and Vilja in As to the Ansatz F(r) = 1 + Ar, 

although we have obtained exact solutions for the field equations, we have found enormous 
difficulty in expressing f(R) in terms of the scalar curvature and cosmological constant for 
d > 4. As an explicit example we have considered d = 4. Moreover, we have also found exact 
solutions for the field equations considering F(r) = F$r and F(r) = For n . For these cases the 
corresponding expression for f(R) have been obtained. A further case that we have also studied 
is for F(r) = 1 + Ar n ; unfortunately this case can only be handled approximately; for this case 
we have given the solutions for u(r) and corresponding scalar curvature R. 

Although the main objective of this paper is the analysis of modified theories of gravity in 
higher dimensions, we also have focused our attention to the case d = 2. Exact solutions for 
different Ansatze have been found. For F(r) = Ar + B, a specific approach, namely considering 
Ar < 1, can avoid logarithmic terms arising in the solution for the metric tensor components. 
These logarithmic terms, on the other hand, could be associated with Newtonian potential which 
we have decided, for simplicity only, do not deal with. For F(r) = For n , with n / 1, we have 
explicitly exhibited exact solutions for u(r), which presents a singular behavior for positive value 
of n. Finally, we have briefly considered the inclusion of the Chern-Simons topological term for 
gravity. In this context we have shown that de Sitter solution is also solution of the complete 
equation of motion. 

Many researcher have studied different f(R) models recently. The purpose of this paper is 
to provide a general approach about how to obtain exact and approximated solutions of the 
field equation considering an arbitrary radial function F(r) in a time independent spherically 
symmetric higher-dimensional spacetime. Although we have consider only vacuum solutions, we 
have reasons to believe that many subjects still deserve to be investigated in higher-dimensional 
spacetime scenario. A natural extension could be to consider the presence of matter fields, to 
derive the field equations and develop a similar analysis as we have done here. 
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A Explicit derivation of the solution for Y{r) =constant 

Here in this Appendix we briefly present some steps adopted in the obtainment of (|26p , solution 
of the field equations (fl~6j) and (JTTJ) for Y(r) = Yq =constant. First, let us focus our attention 
on the corresponding homogeneous equation of (I17p . Defining a new variable z = ^r, we have: 



z 2 {z + l)u"{z)+[z 2 {d 



2) + z(d - 3)] u'(z) - 2 [z(d - 1) + d - 2] u(z) = 



(74) 



whose solution is 




(75) 
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where a and b are integration constants. The above solution is expressed in terms of the Phi- 
function [22j, defined as shown below: 



$(u,s,a) = 



n=0 



(a + n) s ' 



(76) 



which converges for 1 1? | < 1 and a 7^ 0, —1, —2, ... . It is simple to verify the following identity, 
for z > 1: 



1 



l,d - 1 



-d 



i 

V 1 

^ n(-z) n 

n=l 



+ In 1 + 



1 



(77) 



Substituting the above expression into (|75p and rescaling the integration constants, the solution 
can rewritten as 



Uh(z) = az 2 + bz 2 



(78) 



Finally, by recovering the original radial coordinate, r, we can express this solution as follows: 

<i 



Uh(r) = c 2 r 2 + cir^ 



ln ( 1 + £) + S^¥r 



(79) 



The particular solution, u p (r), of the field equations is straightforwardly obtained. So, taking 
into account both solutions, we can write u(r) = u p (r) + Uh(r), which reads 



u(r) 



2 2 

C2f + c\r 



r b_\ (-1)" fB\ n 

11 I + Ar) + ^ nr n \ ~A ) 

x 7 n=l x y 



+ Y 



2 A 2 r 2 ( B\ 
In 1 + — 



d — 1 B 2 



Ar J 



+ 1 



2 Ar' 
d-\~B 



(80) 
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